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Abstract 

Semi-classical theories are approximations to quantum theory that treat some de¬ 
grees of freedom classically and others quantum mechanically. In the usual ap¬ 
proach, the quantum degrees of freedom are described by a wave function which 
evolves according to some Schrodinger equation with a Hamiltonian that depends 
on the classical degrees of freedom. The classical degrees of freedom satisfy clas¬ 
sical equations that depend on the expectation values of quantum operators. In 
this paper, we study an alternative approach based on Bohmian mechanics. In 
this approach the quantum system is not only described by the wave function, but 
with additional variables such as particle positions or fields. By letting the classical 
equations of motion depend on these variables, rather than the quantum expectation 
values, a semi-classical approximation is obtained that is closer to the exact quantum 
results than the usual approach. We discuss the Bohmian semi-classical approxi¬ 
mation in various context, such as non-relativistic quantum mechanics, quantum 
electrodynamics and quantum gravity. The main motivation comes from quantum 
gravity. The quest for a quantum theory is still going on. Therefore a semi-classical 
approach where gravity is treated classically may be an approximation that already 
captures some quantum gravitational aspects. 


1 Introduction 

Quantum gravity is often considered to be the holy grail of theoretical physics. One 
approach is canonical quantum gravity, which concerns the Wheeler-DeWitt equation 
and which is obtained by applying the usual quantization methods (which were so suc¬ 
cessful in the case of high energy physics) to Einstein’s field equations. However, this 
approach suffers from a host of problems, some of technical and some of conceptual 
nature (such as finding solutions to the Wheeler-DeWitt equation, the problem of time, 
...). For this reason one often resorts to a semi-classical approximation where gravity 
is treated classically and matter quantum mechanically [1,2]. The hope is that such an 
approximation is easier to analyze and yet reveals some effects of quantum gravitational 
nature. 

In the usual approach to semi-classical gravity, matter is described by quantum field 
theory on curved space-time. For example, in the case the matter is described by a 


1 



quantized scalar field, the state vector can be considered to be a function '&(</>) on the 
space of fields, which satisfies a particular Schrodinger equation 

= H{(j),g)'i/{<f),t), (1) 

where the Hamiltonian operator H depends on the space-time metric g. This metric 
satishes Einstein’s field equations 

G^,((7) = 87rG(T|V(</.,5)|'h), (2) 

where the source term is given by the expectation value of the energy-momentum tensor 
operator. 

This semi-classical approximation of course has limited validity. For example, it 
will form a good approximation when the matter state approximately corresponds to a 
classical state (i.e., a coherent state), but will fail to be so when the state is a macroscopic 
superposition of such states. Namely, for such a superposition 'h = ('hi -|- 'h2)/\/2, we 
have that ('h|T^i,|'h) ps 2 \T^u\'^ 2 )'j /2, so that the gravitational field 

is affected by two matter sources, one coming from each term in the superposition. 
However, one expects that according to a full theory for quantum gravity, the states 
I'hi) and |'h 2 ) each have their own gravitational held and that the total state is a 
superposition of those. And, indeed. Page and Geilker showed with an experiment of 
this type that this semi-classical theory is not adequate [2,3]. 

Of course, as already noted by Page and Geilker, it could be that this problem is 
not due to fact gravity is treated classically, but due to the choice of the version of 
quantum theory. Namely, Page and Geilker adopted the Many Worlds point of view, 
according to which the wave function never collapses. However, according to standard 
quantum theory the wave function is supposed to collapse upon measurement. Which 
physical processes act as measurements is of course rather vague and this is the source 
of the measurement problem. But it could be that such collapses explain the outcome of 
their experiment. If an explanation of this type is sought, one should consider so-called 
spontaneous collapse theories, where collapses are objective, random processes that do 
not in a fundamental way depend on the notion of measurement. (See [4] and [5] for 
actual proposals combining such a spontaneous collapse approach with respectively (2) 
and its non-relativistic version.) 

In this paper, we consider another approach to quantum mechanics, namely Bohmian 
mechanics [6-9]. Bohmian mechanics solves the measurement problem by introducing an 
actual configuration (particle positions in the non-relativistic domain, particle positions 
or fields in the relativistic domain [10]) that evolves under the influence of the wave 
function. According to this approach, instead of coupling classical gravity to the wave 
function it is natural to couple it to the actual matter configuration. For example, in the 
case of a scalar field there is an actual field cps whose time evolution is determined by 
the wave functional 'k. There is an energy-momentum tensor g) corresponding 

to this scalar field and this tensor can be introduced as the source term in Einstein’s 
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field equations: 


(3) 


GfMi/{g) — , g) ■ 

This approach solves the problem with the macroscopic superposition, since the energy- 
momentum tensor will correspond to just one of the macroscopic matter distributions. 

However, there is an immediate problem with this approach, namely that equation 
(3) is not consistent. The Einstein tensor Gf^i, is identically conserved, i.e., = 0. 

So the Bohmian energy-momentum tensor g) must be conserved as well. How¬ 

ever, the equations of motion for the scalar field do not guarantee this. (Similarly, in the 
Bohmian approach to non-relativistic systems, the energy is generically not conserved.) 

We will explain that the root of the problem seems to be the gauge invariance, which 
in this case is the invariance under spatial diffeomorphisms. Because the scalar field and 
the space-time metric are connected by spatial diffeomorphisms, it seems that one can 
not just assume the metric to be classical without also assuming the scalar field (/)b to 
be classical (in which case the energy-momentum tensor is conserved). 

We will see that a similar problem arises when we consider a Bohmian semi-classical 
approach to scalar electrodynamics, which describes a scalar field interacting with an 
electromagnetic field. In this case, the wave equation for the scalar field is of the form 

= H{(j),A)'^{(j),t), (4) 

where A is the vector potential. There is also a Bohmian scalar field ipB and a charge 
current ^) that could act as the source term in Maxwell’s equations 

d,F^‘'{A)=f{<i>B,A), (5) 

where is the electromagnetic field tensor. In this case, we have dyd^F^^ = 0 due to 
the anti-symmetry of F^’^. As such, the charge current must be conserved. However, the 
Bohmian equation of motion for the scalar field does not imply conservation. Hence, just 
as in the case of gravity, a consistency problem arises. We will find that this problem 
can be overcome by eliminating the gauge invariance, either by assuming some gauge 
or (equivalently) by working with gauge-independent degrees of freedom. In this way, 
we can straightforwardly derive a semi-classical approximation starting from the full 
Bohmian approach to scalar electrodynamics. For example, in the Coulomb gauge, the 
result is that there is an extra current Jq which appears in addition to the usual charge 
current and which depends on the quantum potential, so that Maxwell’s equations read 

d^F^^'^{A)=f{cfB,A)+fQ{f,B,A). ( 6 ) 

While it is easy to eliminate the gauge invariance in the case of electrodynamics, this 
is notoriously difficult in the case of general relativity. One can formulate a Bohmian 
theory for the Wheeler-DeWitt approach to quantum gravity, but the usual formulation 
does not explicitly eliminate the gauge freedom arising from spatial diffeomorphism 
invariance. Our expectation is that one could find a semi-classical approximation given 
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such a formulation. At least we find our expectation confirmed in simplified models, 
called mini-superspace models, where this invariance is eliminated. We will illustrate 
this for the model described by the homogeneous and isotropic Friedman-Lemaitre- 
Robertson-Walker metric and a uniform scalar field. 

In this paper, we are merely concerned with the formulation of Bohmian semi- 
classical approximations. Practical applications will be studied elsewhere. Such appli¬ 
cations have already been studied for non-relativistic systems in the context of quantum 
chemistry [11-16]. It appears that Bohmian semi-classical approximations yield better 
or equivalent results compared to the usual semi-classical approximation. (They are 
better in the sense that they are closer to the exact quantum results.) This provides 
good hope that also in other contexts, such as quantum gravity, the Bohmian approach 
also gives better results. Potential applications might be found in inflation theory, where 
the back-reaction from the quantum fluctuations onto the classical background can be 
studied, or in black hole physics, to study the back-reaction from the Hawking radiation 
onto space-time. 

Other semi-classical approximations have been proposed, see for example [17-20], 
and in particular [21,22] where also Bohmian ideas are used. We will not make a 
comparison with these proposals here. 

The semi-classical approximation is just one practical application of Bohmian me¬ 
chanics. In recent years, others have been explored. See [23, 24] for overviews. For 
example, one may use Bohmian approximation schemes to solve problems in many- 
body systems [25, 26] or one may use Bohmian trajectories as a calculational tool to 
simulate wave function evolution [27]. So even though Bohmian mechanics yields the 
same predictions as standard quantum theory (insofar the latter are unambiguous), it 
leads to new practical tools and ideas. 

The outline of the paper is as follows. We will start with an introduction to Bohmian 
mechanics in section 2. In section 3, we present the Bohmian semi-classical approx¬ 
imation to non-relativistic quantum theory and its derivation from the full Bohmian 
theory. Then we will respectively discuss the Bohmian semi-classical approximation to 
the quantized Abraham model in section 4, scalar quantum electrodynamics in section 
5 and quantum gravity in section 6. In the latter section, we will work out a simple 
concrete example to compare the usual semi-classical approximation to the Bohmian 
one and find that the latter gives better results. 

2 Bohmian mechanics 

Non-relativistic Bohmian mechanics (also called pilot-wave theory or de Broglie-Bohm 
theory) is a theory about point-particles in physical space moving under the influence of 
the wave function [6-9]. The equation of motion for the configuration X = (Xi,..., X„) 
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( 7 ) 


of the particles is given by^ 


where v'^ = (vf ,..., ) with 


= -Im 

rrik 




= —Vfc5 

ruk 


( 8 ) 


and il^ = |V^|e''^. The wave function • • • > ^n) itself satisfies the non- 

relativistic Schrodinger equation 


idt'ip{x,t) 



1 

2mk 


Vl + V{x) 


i){x,t). 


(9) 


For an ensemble of systems all with the same wave function tp, there is a distinguished 
distribution given by IV'p, which is called the quantum equilibrium distribution. This 
distribution is equivariant. That is, it is preserved by the particles dynamics (7) in the 
sense that if the particle distribution is given by |V'(x,to)P at some time to, then it is 
given by \'ip{x,t)\‘^ at all times t. This follows from the fact that any distribution p that 
is transported by the particle motion satisfies the continuity equation 


dtp + '^'^k ■ i^tp) = 0 

k=l 


( 10 ) 


and that IV'P satisfies the same equation, i.e., 


+j;Vfc-(vf|V^|2) = 0, 

fc=i 


( 11 ) 


as a consequence of the Schrodinger equation. It can be shown that for a typical initial 
configuration of the universe, the (empirical) particle distribution for an actual ensemble 
of subsystems within the universe will be given by the quantum equilibrium distribution 
[28, 29]. Therefore for such a configuration Bohmian mechanics reproduces the standard 
quantum predictions. 

Note that the velocity field is of the form where = (jf,...,jn) with 

= lm{'i(j*'Vk'tl’)/mk is the usual quantum current. In other quantum theories, such 
as for example quantum field theories, the velocity can be defined in a similar way by 
dividing the appropriate current by the density. In this way equivariance of the density 
will be ensured. (See [30] for a treatment of arbitrary Hamiltonians.) 

This theory solves the measurement problem. Notions such as measurement or 
observer play no fundamental role. Instead measurement can be treated as any other 
physical process. 

There are two aspects of the theory that are important for deriving the semi-classical 
approximation. Firstly, Bohmian mechanics allows for an unambiguous analysis of the 

^Throughout the paper we assume units in which h = c = 1. 
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classical limit. Namely, the classical limit is obtained whenever the particles (or at least 
the relevant macroscopic variables, such as the center of mass) move classically, i.e., 
satisfy Newton’s equation. By taking the time derivative of (7), we find that 

mfcXfc(t) = -Vk{V{x) + , (12) 

where 

^ 2mfc \^p\ 

is the quantum potential. Hence, if the quantum force —^kQ^ is negligible compared to 
the classical force — then the /c-th particle approximately moves along a classical 
trajectory. 

Another aspect of the theory is that it allows for a simple and natural definition 
for the wave function of a subsystem [8,28]. Namely, consider a system with wave 
function ^(rE,y) where x is the conhguration variable of the subsystem and y is the 
configuration variable of its environment. The actual conhguration is {X,Y), where X 
is the conhguration of the subsystem and Y is the conhguration of the other particles. 
The wave function of the subsystem x{x, t), called the conditional wave function, is then 
dehned as 

X{x,t) = i^{x,Y{t),t). (14) 

This is a natural dehnition since the trajectory X{t) of the subsystem satishes 

Xit) = v^{X{t),Y{t),t) = v^iXit),t). (15) 

That is, for the evolution of the subsystem’s conhguration we can either consider the 
conditional wave function or the total wave function (keeping the initial positions hxed). 
(The conditional wave function is also the wave function that would be found by a 
natural operationalist method for dehning the wave function of a quantum mechanical 
subsystem [31].) The time evolution of the conditional wave function is completely 
determined by the time evolution of ■0 and that of Y. This makes that the conditional 
wave function does not necessarily satisfy a Schrodinger equation, although in many 
cases it does. This wave function collapses according to the usual text book rules when 
an actual measurement is performed. 

We will also consider semi-classical approximations to quantum held theories. More 
specihcally, we will consider bosonic quantum held theories. In Bohmian approaches 
to such theories it is most easy to introduce actual held variables rather than particle 
positions [10, 32]. To illustrate how this works, let us consider the free massless real 
scalar held (for the treatment of other bosonic held theories see [32]). Working in the 
functional Schrodinger picture, the quantum state vector is a wave functional 'I'(0) 
dehned on a space of scalar helds in 3-space and it satishes the functional Schrodinger 
equation 

idt'^{(j),t) = ^ j d^x -F V0(x) • V0(x)^ ^(0T)- (16) 
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The associated continuity equation is 


dt\^{(t>,t)\^ + J = 0 ’ 

where T = |'I'|e‘'^. This suggests the guidance equation 


(17) 




SS{cl>,t) 

5(j){x) 


<A(x)=<i>(x,t) 


(18) 


(Note that in this case we did not notationally distinguish the actual field variable from 
the argument of the wave functional.) Taking the time derivative of this equation results 
in 


where 


□0(x,t) = - 


Q^ = - 




S(p{x) 


2T 


d-^x 


0(x) = (/>(x,t) 

(19) 


(20) 



where Q'^ is the quantum potential. The classical limit is obtained whenever the quan¬ 
tum force, i.e., the right-hand side of equation (19), is negligible. Then the field approx¬ 
imately satisfies the classical field equation □(/> = 0. 

One can also consider the conditional wave functional of a subsystem. A subsystem 
can in this case be regarded as a system confined to a certain region in space. The 
conditional wave functional for the field confined to that region is then obtained from 
the total wave functional by conditioning over the actual field value on the complement 
of that region. However, in this paper we will not consider this kind of conditional wave 
functional. Rather, there will be other degrees of freedom, like for example other fields, 
which will be conditioned over. 

This Bohmian approach is not Lorentz invariant. The guidance equation (18) is 
formulated with respect to a preferred reference frame and as such violates Lorentz 
invariance. This violation does not show up in the statistical predictions given quantum 
equilibrium, since the theory makes the same predictions as standard quantum theory 
which are Lorentz invariant.^ The difficulty in finding a Lorentz invariant theory resides 
in the fact that any adequate formulation of quantum theory must be non-local [33]. One 
approach to make the Bohmian theory Lorentz invariant is by introducing a foliation 
which is determined by the wave function in a covariant way [34]. In this paper, we 
will not attempt to maintain Lorentz invariance. As such, the Bohmian semi-classical 
approximations will not be Lorentz invariant, (very likely) not even concerning the 
statistical predictions. This is in contrast with the usual approach like the one for 
gravity given by (1) and (2) which is Lorentz invariant. However, this does not take 
away the expectation that the Bohmian semi-classical approximation will give better or 
at least equivalent results compared to the usual approach. 


^Actually, this statement needs some qualifications since regulators need to be introduced to make 
the theory and its statistical predictions well defined [32]. 
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3 Non-relativistic quantum mechanics 

3.1 Usual versus Bohmian semi-classical approximation 

Consider a composite system of just two particles. The usual semi-classical approach 
(also called the mean-field approach) goes as follows. Particle 1 is described quantum 
mechanically, by a wave function x(xi,t), which satisfies the Schrodinger equation 

idtx{^i,t)= -^Vl + V{xi,X2{t)) (21) 

2mi 

where the potential is evaluated for the position of the second particle X 2 , which satisfies 
Newton’s equation 

m2X2{t) = -{x ^2C(xi,X2)|^^^X2(t) t) 

= /'(i3xi|x(xi,t)p[-V2P(xi,X2)] . (22) 

J X2=X2d) 

So the force on the right-hand-side is averaged over the quantum particle. 

An alternative semi-classical approach based on Bohmian mechanics was proposed 
independently by Gindensperger et al. [11] and Prezhdo and Brookby [12]. In this 
approach there is also an actual position for particle 1, denoted by Xi, which satisfies 
the equation 

Xi(t)=v>^(Xi(t),t), (23) 

where 

= —Im —— , (24) 

mi X 

and where x satisfies the Schrodinger equation (21). But instead of equation (22), the 
second particle now satisfies 

m2X2(t) = -V2P(Xi(t),X2)|^^^X2(t) ’ 

where the force depends on the position of the first particle. So in this approximation 
the second particle is not acted upon by some average force, but rather by the actual 
particle of the quantum system. This approximation is therefore expected to yield a 
better approach than the usual approach, in the sense that it yields predictions closer to 
those predicted by full quantum theory, especially in the case where the wave function 
evolves into a superposition of non-overlapping packets. This is indeed confirmed by a 
number of studies, as we will discuss below. 

Let us first mention some properties of this approximation and compare them to the 
usual approach. In the mean field approach, the specification of an initial wave function 
x{x,to), an initial position ^ 2 (^ 0 ) and velocity X 2 {to) determines a unique solution for 
the wave function and the trajectory of the classical particle. In the Bohmian approach 
also the initial position Xi (to) of the particle of the quantum system needs to be specified 
in order to uniquely determine a solution. Different initial positions Xi{to) yield different 
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evolutions for the wave function and the classical particle. This is because the evolution 
of each of the variables Xi,X2,x depends on the others. Namely, the evolution of x 
depends on X 2 via (21), whose evolution in turn depends on Xi via (25), whose evolution 
in turn depends on x via (23). (This should be contrasted with the full Bohmian theory, 
where the wave function acts on the particles, but there is no back-reaction from the 
particles onto the wave function.) 

The initial configuration Xi(to) should be considered random with distribution 
|x(x, to)P- However, this does not imply that Xi(t) is random with distribution |x(x, t)p 
for later times t. It is not even clear what the latter statement should mean, since dif¬ 
ferent initial positions Xi(to) lead to different wave function evolution; so which wave 
function should x(x, t) be? 

This semi-classical approximation has been applied to a number of systems. Prezhdo 
and Brookby studied the case of a light particle scattering off a heavy particle [12]. 
They considered the scattering probability over time and found that the Bohmian semi- 
classical approximation was in better agreement with the exact quantum mechanical 
prediction than the usual approximation. The Bohmian semi-classical approximation 
gives probability one for the scattering to have happened after some time, in agreement 
with the exact result, whereas the probability predicted by the usual approach does 
not reach one. The reported reason for the better results is that the wave function 
of the quantum particle evolves into a superposition of non-overlapping packets, which 
yields bad results for the usual approach (since the force on the classical particle contains 
contributions from both packets), but not for the Bohmian approach. These results were 
confirmed and further expanded by Gindensperger et al. [13]. Other examples have been 
considered in [11,14,15]. In those cases, the Bohmian semi-classical approximation gave 
very good agreement with the exact quantum or experimental results. It was always 
either better or comparable to the usual approach. These results give good hope that 
the Bohmian semi-classical approximation will also give better results than the usual 
approximation in other domains such as quantum gravity. 

3.2 Derivation of the Bohmian semi-classical approximation 

The Bohmian semi-classical approach can easily be derived from the full Bohmian the¬ 
ory.^ Consider a system of two particles. In the Bohmian description of this system, we 
have a wave function ' 0 (xi,X 2 ,t) and positions Xi{t)^X 2 {t), which respectively satisfy 

®The derivation is very close to the one followed by Gindensperger et al. [11]. A difference is that they 
also let the wave function of the quantum system depend parametrically on the position of the classical 
particle. This leads to a quantum force term in the equation (25) for particle 2. However, this does not 
seem to lead to a useful set of equations. In particular, they can not be numerically integrated by simply 
specifying the initial wave function and particle positions. In any case, Gindensperger et al. drop this 
quantum force when considering examples [11,13,14], so that the resulting equations correspond to the 
ones presented above. 
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the Schrodinger equation 


+ V'(x.,X 2 )] ^ (26) 

and the guidance equations 

Xi (t) = vf (Xi (t), X 2 (t), t), X 2 (t) = 4 (Xi {t ), X 2 (t), t). (27) 

The conditional wave function x(xi,t) = '0(xi, X 2 (t), t) for particle 1 satishes the equa- 
tion ^ 

\dtx{y^i,t)= + V"(xi,X 2 (t))^ X(xi,t)+/(xi,t), (28) 

where 

^(xi,t) = f- 7 ^V'(xi,X 2 ,t)) +iV 2 V’(xi,X 2 ,t) • V^(Xi(t),X 2 (t),t) . 

^ ^ x.=X.(t) 

(29) 

So in case I is negligible in (28), up to a time-dependent factor times xf" we are led to 
the Schrodinger equation (21). This will for example be the case if m 2 is much larger 
than mi {I is inversely proportional to m 2 ) and if the wave function slowly varies as a 
function of X 2 - We also have that 

m2X2(t) = -V2[F(Xi(t),X2) + g’^(Xi(f),X2,t)] , (30) 

X2=X20) 

with the quantum potential. We obtain the classical equation (25), if the quantum 
force is negligible compared to the classical force. 

In this way we obtain the equations for a semi-classical formulation. In addition, we 
also have the conditions under which they will be valid. For other quantum theories, 
such as quantum gravity, we can follow a similar path to hnd a Bohmian semi-classical 
approximation. 

4 The quantized Abraham model 

The next theory we consider is the Abraham model, which on the classical level de¬ 
scribes extended, rigid, charged particles, which interact with an electromagnetic field 
[35]. We will consider a semi-classical approximation, first considered by Kiessling [36], 
where the charges are treated quantum mechanically, while the electromagnetic field is 
treated classically. Just as in the case of non-relativistic particles, it is straightforward 
to find a consistent set of equations, i.e., no difficulties such as the one mentioned in the 
introduction are met. We will again start from the full Bohmian theory to guide us to 
the semi-classical approximation. 

/ contains a term of the form f{t)x, then it can be eliminated by changing the phase of y by a 
time-dependent term. 
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In the classical Abraham model, the charge distribution of the i-th particle is centered 
around a position qj and is given by /9i(x) = e(/?(x —q^), where the function (/? is assumed 
to be smooth, radial, with compact support and normalized to 1. The particles move 
under the Lorentz force law 


mqi = e [E;p(qi) + q* x B^(qi)] , (31) 

where E,p(x) is shorthand for the convolution (E * (p){x) and similarly for B(^(x). E 
and B are respectively the electric and magnetic field. In terms of the electromagnetic 
potential = (Aq, A), they are given by E = —dtA — VAq and B = V x A. The 
electromagnetic field satisfies Maxwell’s equations 




(32) 


with = d^A'^ — and charge current given by 

n n 

j°(x,t) =/)(x,f) = ^e(/9(x - qi(f)), j(x,f) = ^eqi(t)v9(x-qi(t)). (33) 

i=l i=\ 


Note that although we have employed a covariant notation this model is not covariant, 
due to the rigidity of the particles. 

The quantization of the classical model is straightforward [35]. In the Coulomb gauge 
V • A = 0, the Schrodinger equation reads 


idt^ = 


-T g (V, - ieAj(q.))" + C(,) + / A (-1^ 


+ -(V X A’’(x))- 


|2 ' 2 


(34) 

for the wave function '^{q, A'^), where q = (qi,..., q^), A"^ is the transverse part of the 
vector potential^ and 14 is the Coulomb potential;® 




= -^ j d^xpix)^p{x) 


1 

Svr j 

9 n 




|x-y| 


8. ^ 


v?(x)(/^(y) 


ij=l ' 


x-y + qj-qjl 


(35) 

(36) 

(37) 


So for the electromagnetic field we have used the functional Schrodinger representation 
(see e.g. [32] for more details). 

In the Bohmian approach we have the following guidance equations 

q,(() = 4 [V,S(,, A^, t) - eAj(q.)] . (38) 

®We can write A = + A^, with A'^ = A — V V ■ A and A^ = V - A respectively the 

transverse and longitudinal part of the vector potential. The Coulomb gauge corresponds to A^ = 0. 
®We have used the notation ^ f Fy |^[^y| • 
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A^(x,t) = 


5S{q,A'^,t) 


(5A^(x) |g(i),A7’(x,t) ’ 

where= |'I'|e''^. Taking the time derivative, we find 

mqi = e [E^(qi) + q* x B,^(qi)] - ViQ , 




where 


o = -—— 

2m ^ IV’I 2|V^| 


'd^x^ 

’ 


7- = (o 


are respectively the quantum potential and what can be called the quantum charge 
current, which enters Maxwell’s equation in addition to the usual current. Both currents 
are conserved, i.e., = d^jq = 0. These equations are written in terms of and 

it is assumed that the Coulomb gauge V • A = 0 holds. So the equations (41) are 
equivalent to 

□A^=/+jQ, A^ = 0, V2Ao+j° = 0, VAo=j^, (43) 

where A^ is the longitudinal part of the vector potential. The first equation follows 
from differentiating (39) with respect to time, the second one is implied by the Coulomb 
gauge, and the third (and the fourth) defines Aq in terms of the charge density. 

We can now easily consider a semi-classical limit where either the charges or the 
electromagnetic field approximately behave classically. Let us consider the latter case. 
In that case, we assume the quantum current jg to be negligible in (41). 

Consider further the conditional wave function for the charges xid-, t) = V’(9) A'^(x, t)), 
where (g(t), A^(x, t)) is a particular solution to the guidance equations. It satisfies the 
equation 

V n 

i9tx= + yc{q) X + I, (44) 

i=l 

where 


^ ^^<**>'* + iaA^(x,t). ^^) 

Ad" (x)=Ad’ (x,t) 

(45) 

Whenever I is negligible, up to a time-dependent factor times x (cf. footnote 4), the 
wave equation reduces to 


n 

'^dtX= + 


Together with 


q, = - (V,5 - eAj(q,)) , = f 
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and the Coulomb gauge, this defines a consistent semi-classical approximation. The 
second-order equation for the charges is still of the form (40), but now with quantum 
potential 


Q = -—y 

^ 2m ^ 
1=1 


Ixl 


(48) 


This semi-classical theory was studied by Kiessling [36]. Kiessling also considers 
charges with spin and possible alternative guidance equations. Our analysis shows how 
it may be derived from Bohmian approach to the quantized Abraham model. 


5 Scalar electrodynamics 

So far the semi-classical approximations could have easily been guessed without consid¬ 
ering the full quantum theory. In this section, we consider semi-classical approximations 
to scalar electrodynamics, i.e., a scalar field interacting with the electromagnetic field, 
for which this would have been much harder. The approximation will be obtained by 
starting from different but equivalent formulations of the Bohmian approach to quantum 
electrodynamics. These formulations can be found either by considering different gauges 
or by working with different choices of gauge-independent variables (which more or less 
amounts to the same thing [32,37]). The gauges we will consider here are the Coulomb 
gauge, the unitary gauge and the temporal gauge. However, while the Coulomb and the 
unitary gauge completely fix the gauge freedom, the temporal gauge does not. We will 
find that in order to a (consistent) semi-classical approximation, it will appear necessary 
to completely eliminate the gauge freedom (or at least separate gauge degrees of freedom 
from gauge-independent degrees of freedom). As such, the Bohmian formulation in the 
temporal gauge does not immediately seem to lead to a semi-classical approximation. 

5.1 Coulomb gauge 

5.1.1 Classical theory 

We start by formulating the classical theory of scalar electrodynamics. The Lagrangian 
is 

L = jd^x (^{D^<prDf^cj) - , (49) 

where + ieA^ is the covariant derivative, with = (Aq, A) and = 

d^A’^ — d'^A^. The corresponding field equations are 

+ m^cP = 0, = f , (50) 

where 

f = ie (51) 

is the charge current. The theory has a local gauge symmetry 

</> ^ ^ . (52) 
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In this section we will focus on the Coulomb gauge V-A = 0. Writing A = A^ +A^ 
(cf. footnote 5), this gauge corresponds to A^ = 0 and the equations of motion reduce 
to 

+ 2\eAQdt + ieAo — c^Aq + 2ieA'^ • V + </> = 0 , (53) 

□A^ + Vio=j, (54) 

- V^Ao = JO , (55) 

with j = \e{(j)V(f)* — — 2e^A^|())p and jo = — 2e^AQ|(/)p. Equation 

(55) can be written as 

(V^ - 2e^\^\^) Ao = -ie(</.> - 00*). (56) 

This is not a dynamical equation but rather a constraint on Ao; it could be used to solve 
for Ao in terms of 0. 


5.1.2 Bohmian approach 


The classical theory can easily be quantized using the Coulomb gauge. The same quan¬ 
tum field theory can be obtained by eliminating the gauge degrees of freedom on the 
classical level and then quantizing the remaining gauge-independent degrees of freedom 
[32]. In the resulting Bohmian approach there are actual fields 0 and A^ that are guided 
by a wave functional T(0, A^,t) which satisfies the functional Schrodinger equation^ 


idt^ = 
where 


d^x - 


52 


1^ 1 


5(j)*64> 


+ |(V-ieA^)0|2W|0|"--C^C-- 


1 <52 


1 


2 V2 2 5A^‘^ 2 


+ o(VxA 


T\2 


4', 


C(x) = e ' 0*(x) 


50* (x) 


-0(x) 


50(x) 


(57) 


(58) 


is the charge density operator in the functional Schrodinger picture. The first three 
terms in the Hamiltonian correspond to the Hamiltonian of a scalar field minimally 
coupled to a transverse vector potential. The fourth term corresponds to the Coulomb 
potential and the remaining terms to the Hamiltonian of a free electromagnetic field. 

The guidance equations are 




6S 

Scj)* 



cs, 



(59) 


with T = I'hje*'^. Defining 

Ao = -i^CS, (60) 

^The wave functional should be understood as a functional of the real and imaginary part of 0. In 
addition, writing <j> = (0r + we have that the functional derivatives are given by Wirtinger 

derivatives 5/(50 = {SI5<f>r — i5/50i)/\/2 and 5/50 = (5/50r -I- i5/50i)/\/2. 
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we can rewrite the guidance equation for the scalar field as 


-Do</> 


5S 

5(j)* 


( 61 ) 


With this definition of the classical equation (56) is satisfied. 

Taking the time derivative of the guidance equations we obtain, after some calcula¬ 
tion, that 

D,DK4,-m^4, = -^, (62) 

(63) 

where the Coulomb gauge is assumed. The quantum potential is given by 


Q 



d^x 





1 <52 \ 

26A^) 




(64) 


The current is given by the classical expression (51) and jg = (0, jg), with 


can be considered an additional quantum current. The total current j'^+jq is conserved, 
i.e., + jq) = 0, as a consequence of (62). This is required for consistency of 

(63), since = 0 (due to the anti-symmetry of The current j" satisfies 

= —iCQ and is hence not necessarily conserved. 


5.1.3 Semi-classical approximation: classical electromagnetic field 


Let us first consider a semi-classical approximation by treating the electromagnetic field 
classically and the scalar field quantum mechanically. Similarly as before, it can be found 
by considering the conditional wave functional xi'P} A^(f), t) for the scalar field 

and conditions under which the electromagnetic field approximately behaves classically. 
We will not carry through this procedure, but just give the resulting equations. 

The wave functional xi 4 > 7 1) satisfies the functional Schrodinger equation for a quan¬ 
tized scalar field moving in an external classical transverse vector potential and Coulomb 
potential: 

^dtX = J + |(V - ieA ^)(/>|2 - Km 2|((>|2 - ^C^C^x, ( 66 ) 


where C is defined as before. The actual scalar field satisfies 


Dq(I) 


5S 
6^ ’ 


(67) 


where Aq is defined as before, with S now the phase of x- The vector potential = 
(ylo) A"^) satisfies Maxwell’s equations 




f + fg 


( 68 ) 
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where jg = (0, jg) is an additional quantum current, with 






5(p*S(j) ^ 2‘^V2^) ■ 


So the equations (66)-(68) define the semi-classical approximation. 

We still have that df^{j^^+jq) = 0, so that the Maxwell equations (68) are consistent. 
The correct source term in Maxwell’s equations was found by considering the Bohmian 
approach to the full quantum theory. Without this it would have been hard to guess the 
right current. Just using the current would yield an inconsistent set of equations since 
duj'^ = —iCQ. Adding an extra current to so that the total one would be conserved 
would still leave an ambiguity since one could always add a vector that is conserved. 

Another essential ingredient in our derivation was that the gauge freedom was elim¬ 
inated. We will see in section 5.3 that without such an elimination it does not seem 
possible to derive a semi-classical approximation. 

5.1.4 Semi-classical approximation: classical scalar field 

We can also consider a semi-classical approximation where the scalar held is treated clas¬ 
sically and the electromagnetic held quantum mechanically. In this case, the functional 
Schrodinger equation for the wave function y(A^,t) is 


^dtx = J + X + ieA 


*V</>-</.V<(>*) + e2A^2|^nx. (71) 


The actual potential A^ satishes the guidance equation 


A^ = 


and the scalar held satishes the classical equation 

D^D^cj) — = 0 ) 

where again = (Aq, A^) with Aq dehned by (56). 


5.2 Unitary gauge 

There are other possible semi-classical approximations one can consider. We have 4 
independent held degrees of freedom, namely the two transverse degrees of freedom of 
the vector potential and the two (real) degrees of freedom of the scalar held, and any 
of these or combinations thereof could in principle be assumed classical. Here, we will 
consider two different semi-classical approximations. One where the amplitude of the 
scalar held is assumed classical and another one where all degrees of freedom but the 
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amplitude are assumed classical. While we could consider these approximations using 
the Bohmian formulation in terms of the Coulomb gauge, it is more elegant to consider it 
in a different yet equivalent formulation which can be obtained by imposing the unitary 
gauge. 

5.2.1 Classical theory 

The unitary gauge is given hy (j) = (p*. Writing 

(/) = ?7e'V\/2, (74) 

with rj = \/2|';^| (which can be done where (p 0), this gauge amounts to 0 = 0. The 
classical equations become 

□?? + m^ri — e^A^A^r] = 0 , = j'^ , (75) 

where now = —e^rfA'^. Maxwell’s equations become 

□A + VV • A + e^^A + Vio = 0 , (76) 

(V^ - Ao + V • A = 0. (77) 

The last equation is again a constraint rather than a dynamical equation; it can be used 
to express Aq in terms of A. 




where = (OJq), with 


5A ’ 


(82) 


JQ = - 


and 


Q 


1 




+ 


,52 

5A2 


+ 


g2jy2 





(83) 


This Bohmian approach is equivalent to the approach considered in the previous 
section. This can easily be checked by using the field transformation ((/>, A^) —(r/, A) 
which is obtained by using the polar decomposition (74) for (j) and A = A^ + (with 
inverse transformation 6 = e^V • A).® 


5.2.3 Semi-classical approximation: classical electromagnetic field 


In the case the vector potential A evolves approximately classically, we have the following 
semi-classical approximation. The wave functional xiv^ ^) satisfies 

^dtX = ^ j + (84) 

and the guidance equation 



(85) 


The electromagnetic field satisfies 




( 86 ) 


The appearance of the term containing Aq in the Schrodinger equation requires some 
explanation. As before, the equation can be obtained by assuming certain terms negli¬ 
gible when considering the time derivative of the conditional wave functional xiv^ ^) = 
'I'(r/, A(f), t). In this case, idtx contains the term 


- / d-^x 


2e2r/2(x) 


5S 


(5A(x) 


X = 


I A(x,t) 




-Al{x,t)+ Al{x,t)- 


e^?7^(x) 


SS 


(5A(x) 


A(x,t)^ 


X- (87) 


The term within the big curly bracket on the right hand side would be zero if we 
evaluated it for the actual field r/(x, t) because of (80). If we do not evaluate it then it is 
not necessarily zero. In our semi-classical approximation we assume this term negligible. 
As such, the resulting Schrodinger equation (84) corresponds to the one of a quantized 
scalar field minimally coupled to a classical electromagnetic field. 

Note that there is no consistency issue in this case. Just as in the classical or the full 
Bohmian case, the equation = d^(—= 0 does not follow from the equation 
of ry, but from the Maxwell equations themselves. 

®The equivalence of the Schrodinger picture for the Coulomb gauge and the unitary gauge was 
considered before in [38] (but seems to require a small correction for the kinetic term for ly). 
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5.2.4 Semi-classical approximation: classical scalar field 

In the case we assume the scalar field classical, we have the following semi-classical 
approximation. The wave functional satisfies 

= +(Vx A)2 + eVA2)x (88) 

and guides the field A through the guidance equation 

A = SS/6A - VAo , (89) 

where Aq is defined as in (80). The field rj satisfies the classical equation 

□?7 + m^r] — e^A^A^rj = 0 . (90) 

The Schrodinger equation corresponds to a spin-1 field with mass squared [32]. 

5.3 Scalar electrodynamics: temporal gauge 

In this section, we consider a Bohmian approach to scalar electrodynamics where not 
all gauge freedom is eliminated. It will appear that we need to deal with the remaining 
gauge freedom in order to get an adequate Bohmian semi-classical approximation. 


5.3.1 Classical theory 

The temporal gauge is given by Aq = 0. It does not completely fix the gauge. There is 
still a residual gauge symmetry given by the time-independent transformations : 


e'^V , A ^ A + V0 , 

with 6 = 0. 

In this gauge, the classical equations of motion read 

^ — D^(/) -|- m^cj) = 0 , DA -|- VV • A = j , —V • A = jo , 
where D = V — ieA, j = ie((/>D*(/)* — (jADcj)) and jo = ie((/>*(^ — 


(91) 


(92) 


5.3.2 Bohmian approach 

Quantization of the classical theory leads to the Schrodinger equation® 


id/T = / d^xi — 




+ iDcpl + - 


1 ,52 


1 


2JA^ (93) 


®In [39] this equation is considered to study the semi-classical approximation in the context of stan¬ 
dard quantum theory. 
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so that, assuming the gauge Aq = 0, (92) and (99) can be written as 

= d^F^^=f+fQ, ( 100 ) 

where = (0,jQ). 

This Bohmian approach is equivalent to the one formulated using the Coulomb gauge 
(and hence also to the one formulated using the unitary gauge). To see this, consider 
the field transformation (/>, A —)• defined by (j)' = (/>exp(—ie^V • A) and 

the usual decomposition of A into transverse and longitudinal part. In terms of the 
new variables, the constraint (94) reads 5'^/5A^ = 0, i.e., the wave functional does not 
depend on A^, just on (/> and A^. The Schrodinger equation (93) reduces to the one 
in the Coulomb gauge given in (57). (The latter equivalence is discussed in detail in 
[38].) The guidance equations (95) also reduce to the ones in the Coulomb gauge. They 
yield the extra equation A^ = 0, which implies that the field A^ is static. But A^ is 
the gauge degree of freedom that remains after imposing the temporal gauge. So one 
could gauge-fix it to be zero or just discard it. Therefore, these Bohmian formulations 
are equivalent. (See [32] for a similar comparison in the case of the free electromagnetic 
field.) 
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5.3.3 Usual semi-classical approximation 

In the framework of standard quantum theory, there is a natural semi-classical ap¬ 
proximation that treats the vector potential classically and the scalar field quantum 
mechanically. The scalar field is described by a wave functional t) which satisfies 

'^dtX = J d^x(^ - ( 101 ) 

and the electromagnetic field satisfies the classical Maxwell equations 

= ixlllx) , ( 102 ) 

where 

f / (5T 

= e j 

{x|j|x> = ie(^.D-0* - 0-D^.) , (103) 

and where of course the temporal gauge is assumed. 

This theory is consistent since d^{x\j^\x) = 0, as a consequence of the Schrodinger 
equation (101). It is also invariant under the time-independent gauge transformations 
A ^ A' = A + V0, T((/>) ^ T'((/>) = T(e-‘"V)- 



5.3.4 Bohmian semi-classical approximation 


A natural guess for a Bohmian semi-classical approximation similar to the usual one is 
the following. An actual held 4* is introduced that satishes (f) = 5S/6(p*, where the wave 
functional satishes (101), and the Maxwell equations are = j'^, where is the 

classical expression for the charge current. However, the second-order equation for the 
Bohmian held is 


(f) — D'^cj) -|- rn?4> 


64)* ’ 


(104) 


where Q — J ^ consequence, we have that = —iCQ and hence 

Maxwell’s equations imply that CQ = 0 oi Q = Q{\4)\'^). This is a constraint on the wave 
functional that was absent in the usual semi-classical theory. It also seems to be a rather 
strong condition. It will for example be satished if the scalar held evolves classically 
(i.e., when the right-hand side of (104) is zero) but it is unclear whether there are other 
solutions. 

We arise to a similar conclusion from a more careful approach trying to derive the 
semi-classical approximation from the full Bohmian theory. If we want a semi-classical 
approximation with A classical, then we should require that jg = 0 (since we want 
the wave functional to depend solely on the scalar held and no longer on the vector 
potential). However, due to the constraint (94) (or again using the equation of motion 
(100) for the scalar held), this implies that CQ = 0. 
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So the conclusion seems to be that if we assume A classical, then (p should also 
behave classically. This is not surprising since the gauge symmetry implies that the 
physical (i.e., gauge invariant) degrees of freedom are some combination of the fields A 
and 4>. So one can not just assume A classical and keep (p fully quantum. 

A possible way to develop a semi-classical approximation is thus to separate gauge 
degrees of freedom from gauge-independent ones and assume only some of the latter 
to be classical. One way to do this is as follows. Writing A = A^ -|- we can 
assume A^ to behave classically and not A^, since A^ does not change under a gauge 
transformation, whereas does. We could fully eliminate A^ by introducing the field 
variable cp' = ())exp(—ie^V • A), which would lead to the full Bohmian approach and 
the semi-classical approximation of section 5.1.3. However, we can also formulate a 
semi-classical approximation by keeping the variable A^, but which is still equivalent 
to the one of section 5.1.3. The functional Schrodinger equation is 

i8,4. = y ^ + |D^|2 + (105) 


and the constraint 




= 0 . 


The guidance equations are 


5S ■L_ 

S(j)* ’ 5A^ 


The equation of motion for A^ is 


□A^ = . 


(106) 


(107) 


(108) 


Together these equations imply 

D,D>^P-m^<P = -^, d^F>^'^=r + fQ (109) 

in the temporal gauge, where jg = (0, jg), with 

and quantum potential 

^ ~]¥\ I ^ ^ ^ 2(5A“) 

We still have gauge invariance under time independent gauge transformations. To show 
that this semi-classical approximation is equivalent to the one of section 5.1.3, one just 
has to apply the field transformation {(f>, A^) —((/>', A^) with p' = i;Aexp(—ie^ V • A^) 
(cf. the last paragraph of section 5.1.3). 
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6 Quantum gravity 

6.1 Canonical quantum gravity 

In canonical quantum gravity, the state vector is a functional of a spatial metric /ijj(x) 
on a 3-dimensional manifold and the matter degrees of freedom, say a scalar field 
The wave functional is static and merely satishes the constraints [2]: 

n^{h,^) = 0, ( 112 ) 

n^'^’ih,<P) = 0. (113) 

Their explicit forms are not important here. The latter constraint expresses the fact 
that the wave functional is invariant under inhnitesimal diffeomorphisms of 3-space. 
The former equation is the Wheeler-DeWitt equation. It is believed that this equation 
contains the dynamical content of the theory. However, it is as yet not clear how this 
dynamical content should be extracted. This is the problem of time [2,41]. 

In the Bohmian approach, there is an actual 3-metric and a scalar field, whose dy¬ 
namics depends on the wave function [42-44]. The dynamics expresses how the Bohmian 
configuration changes along a succession of 3-dimensional space-like surfaces.^® Al¬ 
though the wave function is stationary, the Bohmian conhguration will change along 
these surfaces for generic wave functions. This is how the Bohmian approach solves the 
problem of time. 

The structure of the theory is similar to that of scalar electrodynamics in the tempo¬ 
ral gauge, which was discussed in section 5.3. Namely, in both cases there is a constraint 
on the wave functional which expresses invariance under infinitesimal gauge transfor¬ 
mations: spatial diffeomorphisms in the case of gravity and phase transformations in 
the case of scalar electrodynamics. Therefore we may encounter similar complications 
in developing a consistent Bohmian semi-classical approximation. Indeed, as we will 
see below, the Bohmian energy-momentum tensor will not be covariantly conserved and 
hence can not be used in the Einstein equations. This feature is analogous to what we 
saw in the case of scalar electrodynamics. In that case, the Bohmian charge current was 
not conserved so that it could not enter as the source in Maxwell’s equations. The possi¬ 
ble solution to the problem is presumably similar to that in the case of electrodynamics, 
namely the gauge invariance should be eliminated by working with gauge invariant de¬ 
grees of freedom or by choosing a gauge. However, this is a notoriously hard problem in 
the case of gravity. It can be solved in the case of simplified models of quantum gravity 
called mini-superspace models. We will discuss such a model in section 6.4. 

First, we will consider the problem in formulating a Bohmian semi-classical approx¬ 
imation in more detail. Instead of starting from the Bohmian formulation for canonical 

^°The succession of the surfaces is determined by the lapse function and different choices of lapse 
function lead to a different Bohmian dynamics. This is analogous to the fact that in Minkowski space- 
time, the Bohmian dynamics (at least in the usual formulation) depends on the choice of a preferred 
reference frame or foliation. 
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quantum gravity, we will start from a quantized scalar field on an classical curved space- 
time. 


6.2 Semi-classical gravity 


The formulation of quantum field theory on a classical curved space-time in the func¬ 
tional Schrodinger picture was detailed in [45,46]. We assume that the space-time 
manifold A4 is globally hyperbolic so that it can be foliated into space-like hypersur¬ 
faces. A4 is then diffeomorphic to M x S, with S a 3-surface. We choose coordinates 
= (t, x) such that the time coordinate t labels the leaves of the foliation and x are 
coordinates on S. In terms of these coordinates the space-time metric and its inverse 
can be written as 




-Ni -hij) ’ 


9 ^'' = 




(114) 


where N is the lapse function and W = hijN^ are the shift functions, hij is the induced 
Riemannian metric on the leaves of the foliation. The unit vector field normal to the 
leaves is = {N, 0, 0,0). 

For a mass-less scalar field, the Schrodinger equation for this space-time background 
reads 


where T is a functional on the space of fields (/> on S and 




(115) 


(116) 


= + ■ ( 117 ) 

This equation describes the action of the classical metric onto the quantum field. 
The usual way to introduce a back-reaction is by using the expectation value of the 
energy-momentum tensor operator as the source term in Einstein’s field equations, i.e., 

G^,(5) = ('k|V(</>,5)|'k). (118) 


This equation is consistent since V^(T|r^jy((/), (/)|T) = 0. 

In the Bohmian approach, there is an actual scalar field which satisfies the guidance 
equation 

+ (119) 

or, equivalently, 

, 120 ) 

This dynamics depends on the foliation (unlike the wave function dynamics). That is, 
different foliations will lead to different evolutions of the scalar field. 
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Taking the time derivative of this equation, we obtain 


= 


1 5Q 
\/^ ’ 


( 121 ) 


with 


Q = - 


2 m 


(VX 


N (52|T| 

Wfi^\ 


( 122 ) 


the quantum potential. 

To find the energy-momentum tensor for the Bohmian field, we consider the La- 
grangian 

1 r . 

(123) 




from which the equation of motion (121) can be derived. We can then use the usual 
definition of the energy-momentum tensor 


2 5L 


(124) 


to obtain 

. (125) 

This seems to be the natural expression for the Bohmian energy-momentum tensor, 
which could be used in the classical Einstein field equations.However, by taking the 
covariant derivative, we find 


V'%. = -^(^M+2V-^). (120) 

which is generically different from zero, so that Einstein’s equations would not be con¬ 
sistent with this energy-momentum tensor as matter source. 


6.3 Bohmian analogue of the Schrodinger-Newton equation 


In the non-relativistic domain, one often describes the coupling between quantum matter 
and classical Newtonian gravity by the Schrodinger-Newton equation (see e.g. [48]). Eor 
one particle, this equation, which can be derived from the equations (115) and (118) 
[49], is given by 


i(9tV'(x,t) 


-F m<I>(x, t) 
2m 


V’(x,t), 


(127) 


where is the gravitational potential, which satisfies 


V24> = 47rGmlV^l2. (128) 

a more precise derivation, we should start from the Bohmian theory of quantum gravity, deduce 
the modified Einstein field equations (along the lines of [47]), and find the energy-momentum tensor for 
the matter field. 


25 



So the mass distribution that generates the gravitational potential is given by 
This leads to the non-linear Schrodinger equation 

i5tV’(x,t)= f V'(x,t). (129) 

[2m J |x-y| J 

One could consider a Bohmian analogue of the Schrodinger-Newton equation, where 
in addition to (127) we have the guidance equation 

X(t) = v^(X(t),t) (130) 

and the Poisson equation 

V2$(x,t) = 47rGm(5(x-X(t)). (131) 

So in this case the gravitational potential is generated by the Bohmian particle rather 
than the wave function. This leads to the modified Schrodinger equation 

The generalization to many particles is straightforward. In this case, the gravita¬ 
tional potential generated by the Bohmian particles is 

so that the potential in the many-particle Schrodinger equation is given by 

V{x) = '^ mfc4>(xfc, t). (134) 

k 

Such a Bohmian version of the Schrodinger-Newton equation does not immediately 
seem to follow from a more fundamental Bohmian theory. Of course, to possibly obtain 
such an equation, we should introduce positions rather than a field configuration for the 
matter, which could be naturally done for the Dirac field. (It seems that Bohmian field 
ontologies are well-suited for bosonic field theories, while particle ontologies are well- 
suited for fermionic field theories [10]). However, the reason why we think (132) does 
not follow from such a theory is that we obtained nothing of this sort for the Bohmian 
semi-classical Abraham model considered in section 4. In the classical Abraham model, 
the Poisson equation expresses the electric potential in terms of the charges. This may 
suggest a Bohmian analogue in which the electric potential is generated by Bohmian 
point charges. However, this is not what we find in the Bohmian semi-classical approach, 
cf. (46). This being said, it might still be worth investigating the Bohmian analogue of 
the Schrodinger-Newton equation, perhaps as a possible alternative to collapse models. 
(A model similar in spirit was studied in [50].) 
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6.4 Mini-superspace model 

In this section, we consider a symmetry-reduced model of quantum gravity where homo¬ 
geneity and isotropy are assumed. In this model, the spatial diffeomorphism invariance 
is eliminated and we can straightforwardly develop a Bohmian semi-classical approxi¬ 
mation. 

In the classical mini-superspace model, the universe is described by the Friedmann- 
Lemaitre-Robertson-Walker (FLRW) metric 

ds^ = — a(t)^dH|, (135) 

where N is the lapse function, a = e“ is the scale factor^^ and dill is the metric on 
3-space with constant curvature k. Assuming matter that is described by a scalar field 
(/>, the Lagrangian is [51,52] 

where k = S/AttG. 

VG = -he-^^ + ^A (137) 

is the gravitational potential, with A the cosmological constant, and Vm is the potential 
for the matter field. The corresponding equations of motion are, after imposing the 
gauge N = id^ 

+ Vg , (138) 

^ 3q!(^ -|- — 0 . (139) 

(The second-order equation for a which arises from variation with respect to a is re¬ 
dundant since it can be derived from the other two equations.) 

Using the canonical momenta 

TTat = 0 , TTa = , (140) 

we can pass to the Hamiltonian formulation. This leads to the Hamiltonian constraint 
(which is just eq. (138)) 

^^The reason for introducing the variable a is that it is unbounded, unlike the scale factor, which 
satisfies a ^ 0. 

^^The theory is time-reparamaterization invariant. Solutions that differ only by a time- 
reparameterization are considered physically equivalent. Choosing the gauge N = 1 corresponds to 
a particular time-parameterization. 
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Quantization yields the Wheeler-DeWitt equation: 


{Hg + Hm)^^ = 0 , 


(142) 


where 

In the corresponding Bohmian approach [52], there is an actual FLRW metric of 
the form (135) and scalar field, whose time evolutions are determined by the guidance 
equations 

« = = (144) 

where N is an arbitrary lapse function.In the gauge = 1, these equations imply 


-a^ — - ^2*^^ 4- + he + Q'g ^ (145) 

(j) + 3a(j) + d^iVu + + i^Q^) = 0 , (146) 


where 

^ 1 ggldl ^_1 W 

^ 2K^e®" IV’I ’ ^ 2e®“ \ip\ 

We will now look for a semi-classical approximation where the scale factor behaves 
approximately classical. In order to do so, we assume again the gauge N = 1 and we 
consider the conditional wave function = d('/’)«(^))) given a set of trajectories 

Using 

dtx{4>,t) = 9aV’(<(',a)L=„(t)«(^) > (148) 

we can write 

idtx = HmX + I, (149) 

where^® 



I = ^idtx ( d + 
a 


1 




3a^<^^\a{t) J ^ 2Ke^°‘ 


+ - 


^ [id^Sf + iS^^] x+ne^^iVc+Qt) 


a{t) 


a=a[t) 

(150) 

When I is negligible (up to a real time-dependent function times x); (149) becomes the 
Schrodinger equation for a homogeneous matter held in an external FLRW metric. We 


^^Just as the classical theory, the Bohmian approach is time-reparameterization invariant. This is a 
special feature of mini-superspace models [53, 54]. As mentioned before, for the usual formulation of the 
Bohmian dynamics for the Wheeler-DeWitt theory of quantum gravity, a particular space-like foliation 
of space-time or, equivalently, a particular choice of “initial” space-like hypersurface and lapse function, 
needs to be introduced. Different foliations (or lapse functions) yield different Bohmian theories. 

^®To obtain this equation, note that = [(das'll -I- id^S -I- d^lV’l/IV’l]'*/’ + “iidaSdaip, so that 

dai’la^ait) = [(das'll -|- id^ + d^j■i/>|/ji/)|]jc<=c«(t)X + 2idaSdtx/f^- Using this equation together with 
(142) we obtain (149). 
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can further assume the quantum potential Qq to be negligible compared to other terms 
in eq. (145). As such, we are led to the semi-classical theory: 


idtx = HmX , 


1 


dtS + Vg . 


(151) 

(152) 

(153) 


Let us now consider when the term I will be negligible. The quantity in brackets 
in the first term would be zero when evaluated for the actual trajectory (j)(t) (because 
of the guidance equation for a). As such, the first term will be negligible if the actual 
scale factor evolves approximately independently of the scalar field. The second term 
will be negligible if S varies slowly with respect to a or if the term in square brackets 
is approximately independent of (p. In the latter case, the second term becomes a 
time-dependent function times y, which can be eliminated by changing the phase of 
y. Similarly, if Qq <C Vq then the third term also becomes a time-dependent function 
times y. 

In the usual semi-classical approximation, one has (151) and 

=-^{x\Hm\x)+ Vg, (154) 

with y normalized to one. These equations follow from (1) and (2). In the next section, 
we will compare this approximation with the Bohmian one for a particular example. It 
will appear that the latter gives better results than the usual approximation. (Note that 
Vink himself, in his seminal paper [52] on applying the Bohmian approach to quantum 
gravity, considers a derivation of the usual semi-classical approximation, rather than the 
Bohmian one. But he hinted on the Bohmian semi-classical approximation in [55].) 


6.5 Example in mini-superspace 

In this section, we will work out a simple example to compare the Bohmian and usual 
semi-classical approximations to the full Bohmian result. 

We put k: = 1 and assume Vg = Vm = 0. It will also be useful to introduce the time 
parameter r, defined by dTe^°‘ = dt (which corresponds to choosing N = e^“ instead 
of N = 1). Derivatives with respect to r will be denoted by primes. In this way, the 
classical equations (138) and (139) reduce to 

a'2 ^ ^ (j)" = {). (155) 

The possible solutions are 

a = cir-|-C2, 4) = ±cit + C2,, (156) 

where Cj, i = 1,2,3, are constants. In the case ci = 0, the scale factor is constant and 
we have Minkowski space-time. If ci > 0 the universe starts from a big bang and keeps 
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expanding forever. If ci < 0 the universe contracts until a big crunch. If ci 7 ^ 0, the 
corresponding paths in {(j), a)-space are given by 

a = ztcj) + c, (157) 


with c constant. 


6.5.1 Full Bohmian analysis 

In the full quantum case, we have the Wheeler-DeWitt equation 

{dl - dl)i) = 0 

and guidance equations 


a' = -daS, (j)' = dsS 


For the state 


= exp 


m{(j) — a) — 



(158) 

d^S. 

(159) 

{(j) - af 
4(t2 

, (160) 


the guidance equations read a' = 4>' = u, so that we have only classical solutions 

a = UT + Cl , (j) = UT + C2 (161) 


or 


a = 0 + c. 


(162) 


See fig. 1 for some trajectories. 
Similarly, for the state 

V’L (</>,«) = exp 

the solutions are also classical: 

a = VT + Cl 


—iv{(f) + a) — 


i> + af 


= -VT + C2 , 


(163) 


(164) 


or 

a = —(j) + c, (165) 

see fig. 2 . 

Consider now the superposition V’ = i’R + ipL and assume v > u ^ 0. For a —)■ ±00 
the wave functions tpR and ipL are non-overlapping functions of (p so that, asymptotically, 
the Bohmian dynamics is either determined by V’J? or ipR. This means that asymptoti¬ 
cally, i.e., for a —)■ ± 00 , we have classical motion, given either by (162) or (165). Some 
trajectories are plotted in figs. 3 and 4. (Trajectories for the case u = —v were plotted 
in [56].) We see that trajectories starting from the “left”, i.e. (/> < 0, for a —)■ —00 will 
end up on the left, while trajectories that start from the “right”, i.e. </> > 0 , might either 
end up moving to the left or to the right. (If u = v, then trajectories on starting on 
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Figure 1: Some trajectories for a = Figure 2: Some trajectories for 'E'^, cr = 

1 . 1 . 

the left stay on the left and trajectories starting on the right stay on the right.) Some 
trajectories are closed. They correspond to cyclic universes (which oscillate between a 
minimum and maximum scale factor), see fig. 4. 

For trajectories with classical asymptotic behavior, there is possible non-classical 
behavior in the region of overlap, where a ~ 0. For trajectories starting on the left 
there is a transition from a = ucj) to a = —v(f) (which is impossible classically). For some 
trajectories starting on the right, there is the opposite transition. 

Note that there is no natnral measure on the set of trajectories [57], so we can not 
make any probabilistic statements like about the probability for a trajectory to move 
from left to right. 

6.5.2 Usual semi-classical approximation 

Let us first consider the usual semi-classical approximation and compare it to the full 
Bohmian approach: 

idrX = -^dlx, (166) 

= 2{x\Hm\x) =-{x\dl\x) ■ (167) 

The wave equation corresponds to that of a single particle of unit mass in one dimension. 
Hence, {x\d‘^\x) is time-independent and we have that a' is constant. 

The initial wave function which we will use in the semi-classical approximation is 
given by the conditional wave fnnction x(<?i, tq) = Nil>{(p,a{To)) at some time tq, with 
N a normalization constant. The conditional wave functions corresponding to ipR and 


31 





Figure 3: Some trajectories (or pieces 
thereof) for 'I'/j + u = 1, u = 5, 
fj = 1. Three trajectories are high¬ 
lighted which illustrate the possible be¬ 
haviors for trajectories which asymptot¬ 
ically go like a = ±()). 


0.3 

0.2 

0.1 

a 0.0 

- 0.1 

- 0.2 

- 0.3 


Figure 4: Some trajectories (or pieces 
thereof) for u = 1, v = 5, a = 

1. Three trajectories are highlighted. 
Two of these correspond to closed loops, 
which correspond to cyclic universes. 



<P 


ipL are 

Xr(</>,t-o) = NR'il;R{(j),ao) 
Xl{4>,To) = NL'lpL{(t),0!o) 


(27rcr^) ^/“^exp 
(2707^)“^/^ exp 


m{(j) — ao) - 
-w{(j) + ao) 


{(j) - ap)^' 
dcr^ 

{(l) + aof 
4cj2 


(168) 


where ao = a(ro). XR is a Gaussian packet centered around ao and average momentum 
u, while XL is a Gaussian packet centered around —ao and average momentum v. The 
solutions to the Schrodinger equation (166), with these conditional wave functions as 
initial conditions, are given by [7]: 


XR{(t>,r) = [2775^(1-)] ^/“^exp [m((/>R - ur/2) - {(j)R - ur)^/4s(r)cr] , 

Xl{4>, l) = [27rs^(r)]“^/^ exp [-\v{cI)l + ur/2) - (</>l -h ur)^/4s(f )(t] , (169) 

where 


4>R = 4> - cto , </>L = </> + ao , r = r - ro , s(r) = cj( 1 -h ir/2(T^). (170) 


In the following, we assume that 

• ro is small enough, so that ao < 0 and |ao| ^ cr, 

• a'(ro) > 0 ((167) determines a' only up to a sign), 

• u > u > 0 , 
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• 1 ^ au. 


The expectation value (xI-^m|x) is time-independent. So we can calculate it at time 
tq. We have that 

“^{xrIHmIxr) = + ^ ■ ( 171 ) 

Hence a' = \Jv?‘ + Since 1 <C fra, we have that a' k, u, so that we approximately 
obtain the classical solution for the scale factor. 

Similarly, for XL, we have that 

“^{xlIHmIxl) = v'^ + ^ (172) 

and since 1 au < av, we have a' ^ v. 

Now consider the superposition x = {XR + Xl)/V 2- Since XR and xl have ap¬ 
proximately negligible overlap initially, because |ao| a, this state is approximately 
normalized to one. For the same reason we have that 

2(xI^m|x) ~ {xr\Hm\xr) + {xl\Hm\xl) = ^(“^ + ^^) + ^ • (^'^^) 

Hence, for 1 <C au, we have that a' k, {ifi -|- v‘^)l2. As such, the semi-classical 
approximation is very close to the exact result, given that u ~ n. But if u is very 
different from v, the semi-classical approximation is not so good. In particular, we do 
not get the asymptotic behavior that a' = uoi a' = v for early or late times (i.e., r —)■ tq 
or r —)■ oo). 

6.5.3 Bohmian semi-classical approximation 

The equations of motion in the Bohmian semi-classical approximation are 


i5rX = -^9|x, (174) 

4>' = d4>S, (175) 

a'2 = ^'2 _ = _2Q^s . (176) 

Ixl 

For the state V’J?, the solutions of the guidance equation are [7]: 

(p - ao = u{t - To) + {(t)o - ao )^-^—Aflli ^ (I77) 

(7 

where (po = P{ro)- Or using the notation (170): 

pR = uf + , (178) 

a 


where (pR^ = (Pr{to) = (po — ao- With the same assumptions about the constants as in 
the previous section and taking ^ <7, i.e., that the initial value pR^ does not lie 
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too far outside the bulk of the Gaussian packet, we have that \(I)r^q\t/ 2a‘^ <C ut. Using 
in addition that + r^/dcr^ ^ 1 + ^ we hnd that 

(t)R‘^ uf + (t)Rfif/2a^ . (179) 


So we approximately get classical motion and hence the full Bohmian result. 
The classical equation for the scale factor becomes 


/2 2 I 

Oi — Vj -\- 


' Q I / —\ I ' 


r 2 


- 1 


We have that 
3 


(t>R,0UT (pRfi 

' ^ Q I / _\ I “U 


2|5(r)p 2(T^|5(r)| 4(j^ \4(j^|5(r)p 

u\4>R,of\ , 4>\q 


2|s(r)P 2(T3|s(r)| 4(T^ \4cj2|s(r)P 


r 


- 1 


< 


+ 


2|s(r)P 2cj^|s(r)| 4(T^ 


^ 3 u\(t)Rfl\ (pRfi 
- 2u2 + ^2 + 4^4 ’ 


(180) 


r 2 


4iT2|s(r)|2 

(181) 


- 1 


where we used |r|/2iT|s(r)| ^ 1 for the last two terms in order to obtain the last in¬ 
equality. Using the assumptions that 1 <C art and \4>Rfl\ ^ we hnd that 


'2 ^ 2 
a ss u 


(182) 


So similarly as in the case of the usual semi-classical approximation (just with the extra 
condition on the initial value (pRfl), we obtain the full quantum results. For the wave 
function ipR similar results hold. 

Consider now the superposition x = {XR + Xl)/ This superposition corresponds 
to two Gaussian packets that move across each other. Initially and hnally they are 
approximately non-overlapping.^® This means that before and after the wave packets 
cross, the motions of (j) and a are approximately classically, in agreement with the full 
Bohmian results, since they will be determined by either XR oi' XL- This is unlike the 
usual semi-classical approximation. 

A trajectory for the scalar field starting from the left (i.e., (j) <Q) will end up on the 
left, while a trajectory starting on the right will end up on the right. The reason is that 
because of equivariance the probability to start from the left equals the probability to 
end up on the left (which equals 1/2). Since trajectories do not cross (since the dynamics 
is given by a first-order differential equation in time), the probability for trajectories to 
start on the left (right) and end up on the right (left) must be zero. This means that 
for all trajectories starting on the left, we have a transition from a = ucj) to a = —vcj) 
and the opposite transition for trajectories starting on the right. This is unlike the full 
Bohmian analysis where trajectories exist where such a transition does not occur. 

^®The latter statement follows from the fact that the spread of the wave function, which equals 
|s| = aVl + U/4<7h goes like t jla for large time, which is much smaller than the distance between the 
centers of the Gaussians, which equals {u -I- v)f. 
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In conclusion, we see that the Bohmian semi-classical approximation is in better 
agreement with the exact Bohmian results than the usual semi-classical approximation. 
We came to this conclusion by making the comparison on the level of the actual tra¬ 
jectories for a and (j). We did not attempt to make a comparison in the context of 
standard quantum theory, since the standard quantum interpretation of the Wheeler- 
DeWitt equation is problematic due to the problem of time. But it is clear that for 
approaches to quantum theory that would associate (approximately) the classical evolu¬ 
tions (162) and (165) to the superposition (like perhaps the consistent histories 

or many worlds approach) the usual semi-classical approximation would fare worse than 
the Bohmian semi-classical approximation. 

7 Conclusion 

We have shown how semi-classical approximations can be developed using Bohmian 
mechanics. We have obtained these approximations from the full Bohmian theory by 
assuming certain degrees of freedom to evolve approximately classically. This was illus¬ 
trated for non-relativistic systems. If there is a gauge symmetry, like in electrodynamics 
or gravity, then extra care is required in order to obtain a consistent semi-classical 
theory. By eliminating the gauge symmetry (either by imposing a gauge or by work¬ 
ing with gauge-independent degrees of freedom), we were able to find a semi-classical 
approximation in the case of scalar quantum electrodynamics. For quantum gravity, 
eliminating the gauge symmetry (more precisely the spatial diffeomorphism invariance) 
is notoriously hard. We have only considered the simplified mini-superspace approach 
to quantum gravity, which describes an isotropic and homogeneous universe, and where 
the diffeomorphism invariance is explicitly eliminated. More general cases in quantum 
gravity still need to be studied. For example, for the case of inflation theory, where one 
usually considers quantum fluctuations on a classical isotropic and homogeneous uni¬ 
verse, it should not be too difficult to develop a Bohmian semi-classical approximation. 

Apart from possible applications in quantum cosmology, such as inflation theory, it 
might also be interesting to consider potential applications in quantum electrodynamics 
or quantum optics. In particular, since the results may be compared to the predictions 
of full quantum theory, this could give us a handle on where to expect better results 
for the Bohmian semi-classical approximation compared to the usual one in the case of 
quantum gravity where the full quantum theory is not known. That is, it might give us 
better insight in which effects are truly quantum and which effect are merely artifacts 
of the approximation. 

Further developments may include higher order corrections to the semi-classical ap¬ 
proximation. One way of doing this might be by following the ideas presented in [58, 59]. 
As explained there, one might introduce extra wave functions for a subsystem in addi¬ 
tion to the conditional wave function. These wave functions interact with each other and 
the Bohmian configuration. By including more of those wave functions one presumably 
obtains better approximations to the full quantum result. 
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Finally, although we regard the Bohmian semi-classical approximation for quantum 
gravity as an approximation to some deeper quantum theory for gravity, one could also 
entertain the possibility that it is a fundamental theory on its own. At least, there is 
presumably as yet no experimental evidence against it. 
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